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I .  INTRODUCTION 


Over  the  past  several  years,  a  number  of  computer  models  of  the  regenerative 
interior  ballistic  process  have  been  developed.  The  two  models  which  have  been 
used  most  frequently  at  the  Ballistic  Research  Laboratory  are  the  models  developed 
by  Coffee1  and  Gough.2  The  Coffee  model  utilizes  lumped  parameter  descriptions 
for  the  three  main  regions  in  the  regenerative  gun:  the  propellant  reservoir, 
the  combustion  chamber  and  the  barrel.  In  contrast,  Gough  has  chosen  to  treat 
the  reservoir  and  chamber  as  lumped  parameter  regions,  while  utilizing  a 
one -dimensional  flow  model  for  the  barrel  region. 

In  order  to  model  the  barrel  as  a  lumped  parameter  region,  a  pressure 
gradient  model,  analogous  to  the  standard  Lagrange  gradient  used  in  many  solid 
propellant  gun  models,3  is  required.  The  pressure  gradient  model  simulates  the 
drop  in  pressure  from  the  breech  or  barrel  entrance  to  the  projectile  base,  which 
results  from  the  velocity  gradient  developed  as  the  projectile  accelerates  down 
the  barrel.  The  development  of  such  a  pressure  gradient  model  requires  that  the 
unique  features  of  the  regenerative  interior  ballistic  (IB)  process  be  Incorporated 
into  the  model.  The  pertinent  features  are  the  non-zero  gas  velocity  at  the 
entrance  to  the  barrel  and  the  existence  of  a  rarefaction  wave  traveling  along 
the  barrel  toward  the  projectile  base  after  all-burnt. 

In  the  following  paper,  the  development  of  a  modified  Lagrange  gradient 
model,  which  accounts  for  the  unique  features  of  the  regenerative  IB  process, 
is  described.  The  modified  Lagrange  gradient  has  been  incorporated  into  the 
lumped  parameter  interior  ballistic  model  developed  by  Coffee,1  and  computer 
simulations  of  a  hypothetical  120-mm  tank  cannon  have  been  conducted.  Simulations 
were  also  conducted  using  the  Coffee1  model  with  a  standard  Lagrange  gradient 
and  the  Gough2  model  with  a  one -dimensional  barrel  flow  model.  Using  the  results 
from  the  Gough  model  as  a  baseline,  comparisons  of  muzzle  velocities,  pressure 
versus  time  curves,  and  detailed  velocity  and  pressure  profiles  in  the  barrel 
are  presented  and  discussed. 
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II.  THE  LAGRANGE  PRESSURE  GRADIENT 


The  equations  of  motion  governing  the  motion  of  the  gas  in  a  gun  in  the 
region  from  breech  to  projectile  base  are 


ap  ap 

- +  V - 

af  ax 


(l) 


dll  dv 
—  *  V  — 

at  dx 


ap 

'  dx 


with  the  boundary  conditions 


X  BREECH  “  0  • 


(2) 


(3) 


*BA3E-y-  (4) 

^BREECH  *  VL  “  0  ’  (S) 


•W*"*  (6) 
where  u^,  is  the  velocity  of  the  projectile. 

In  the  development  of  the  Lagrange  pressure  gradient  "it  is  assumed  all  the 
propellant  charge,  C,  is  in  gaseous  form  at  the  time  considered".  However,  "the 
theory  applies  without  alteration  if  it  is  assumed  that" ,  prior  to  consumption 
of  all  the  propellant  charge,  "the  unbumt  charge  moves  with  the  gas,  the 
distribution  of  the  solid  along  the  bore  being  the  same  as  the  distribution  of 
gas. "3 


We  now  assume  that  the  density  of  the  gas  (or  gas  plus  unburnt  charge)  is 
uniform  over  the  region  behind  the  projectile, 


2 


We  then  obtain  from  Equation  (1) 


dp  dv 
dt  *dx 


0 


or 


dv  l  dp 
dx~~pJt  ' 

Assuming  a  constant  bore  area,  AB,  (i.e.,  no  chambrage)  and  noting  that  p  -C/A„y 
Equation  (8)  becomes 


dv  1 dy  Up 
dx  y  dt  y 

Integrating  (9)  over  the  region  [0,y],  we  obtain 


(9) 


(10) 


Comer  notes  that  the  term  "Lagrange  approximation"  is  applied  to  Equation  (7) 
or  (10)  and  that  Equation  (7)  leads  to  Equation  (10),  but  it  is  not  true  that 
Equation  (10)  necessarily  implies  Equation  (7). 3 


Substituting  Equation  (10)  in  Equation  (2),  we  have 


1 dP 
p  dx 


or 


dP  Cupf  x_  \ 


(11) 


Integrating  on  [0,y],  and  noting  that 


PCyWBASC 

where  is  the  bore  resistance  pressure  and  Mp  is  the  projectile  mass,  we 


/’(*)-PBasE+2M  ^BASE~^Rts^ 


(-3  • 


For  x  -  0, 


’bIUICH  “  ^BASt  "*■  2M  (  ^  BASS  ~  ^  R£S  )  • 


The  space  mean  pressure  is  defined  by 


-  i  r  v 

P  •  -  I  P{x)dx , 
yJa 


and  upon  substituting  Equation  (12),  we  obtain 


P  “  ^BASt  +  (^BASE  ~  ^*RCs)  * 


The  kinetic  energy  associated  with  the  motion  of  the  gas  is 


KE, 


r  i 

'I,  21-p"' 


dx . 


obtain 


(12) 


(13) 


(14) 


(15) 


(16) 
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Using  Equation  (10)  we  obtain 


KE 


CAS  * 


y 

o 


or 


KECAS-gCu2p  .  (17) 

The  solution  after  "all-burnt"  (to  muzzle  exit)  is  based  on  the  assumption  that 
the  gases  expand  adiabatically  such  that 


-CONSTANT  . 

ALL-BURNT 


However,  this  portion  of  the  solution  is  not  pertinent  to  the  model  developed 
in  subsequent  sections,  and  will  not  be  discussed  here. 


III.  DESCRIPTION  OF  THE  RLPG  INTERIOR  BALLISTIC  PROCESS 


The  RLPG  IB  process  is  based  on  the  injection  of  the  propellant  into  the 
combustion  chamber  during  the  IB  cycle.  The  stages  of  the  IB  process  are  depicted 
in  Figure  1  for  the  RLPG  configuration  known  as  Concept  VI. 

The  system  consists  of  3  regions:  (1)  The  liquid  reservoir  bounded  by  the 
regenerative  piston,  the  fixed  bolt,  and  the  transducer  block;  (2)  the  combustion 
chamber;  and  (3)  the  barrel.  In  regenerative  gun  concepts  similar  to  Concept 
VI,  there  is  a  large  area  change  (about  a  factor  of  4)  from  the  combustion  chamber 
into  the  barrel . 
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Initially,  the  piston  is  seated  on  the  front  of  the  bolt,  sealing  the 
propellant  in  the  reservoir;  see  Figure  1(a).  The  process  is  initiated  by  firing 
an  igniter  into  the  combustion  chamber,  which  generates  an  initial  pressure  in 
the  chamber;  see  Figure  2. 

The  increasing  chamber  pressure  acts  on  the  regenerative  piston,  forcing 
it  to  the  rear,  thus  pressurizing  the  propellant  reservoir.  Due  to  the  area 
difference  across  the  piston  face  from  the  chamber  to  the  reservoir,  the  piston 
acts  as  a  pump,  resulting  in  (1)  a  pressure  in  the  liquid  reservoir  which  is 
higher  than  that  in  the  combustion  chamber,  and  (2)  injection  of  propellant  into 
the  combustion  chamber. 

As  the  piston  moves  to  the  rear  over  the  contoured,  fixed  bolt,  the  injection 
area  first  increases  to  a  maximum  value  associated  with  the  maximum  design 
pressure  of  the  system  and  then  decreases,  decelerating  the  piston  as  it  completes 
its  stroke;  see  Figure  1(b).  When  shot-start  pressure  is  exceeded,  the  projectile 
begins  accelerating  down  the  tube.  Analyses*  of  experimental  data  indicate  that 
early  in  the  ballistic  cycle,  significant  amounts  of  unburnt  propellant  may 
accumulate  in  the  combustion  chamber,  that  this  accumulated  propellant  is  rapidly 
consumed  as  the  pressure  rises  to  its  maximum  value,  and  that  propellant 
subsequently  injected  is  rapidly  consumed  leading  to  very  low  values  of  propellant 
accumulation  in  the  latter  stages  of  the  ballistic  cycle;  see  Figure  2.  Such 
analyses  would  also  indicate  that  over  most  of  the  projectile  travel,  propellant 
is  consumed  rapidly  in  the  combustion  chamber  with  little  accumulation,  and, 
therefore,  only  small  amounts  of  unbumt  propellant  would  be  transferred  into 
the  barrel. 

Utilizing  such  analyses,  most  regenerative  IB  models  have  been  based  on  the 
following  assumptions: 
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©  IGNITER  FUNCTION 
(D  INJECTION  -  LP  ACCUMULATION 
(D  COMBUSTION  OF  LP  ACCUMULATION 
©  QUASI-STABLE  EQUILIBRIUM 


Figure  2.  A  Pressure  Versus  Time  Curve  From  an  Early  Simple 
In-Line  Piston  RLPG  Concept.  Annotated  to 
Show  the  Phases  of  the  Regenerative 
Interior  Ballistic  Process. 

1.  The  combustion  chamber  is  a  homogeneous  region  containing  either  a 
two-phase  mixture  of  propellant  droplets  and  combustion  gases  (if  a  finite  rate 
droplet  burning  model  is  used),  or  combustion  gases  only  (if  an  instantaneous 
propellant  burning  model  is  assumed) . 

2.  Only  combustion  gases  enter  the  barrel  region. 

3.  As  the  gases  enter  the  barrel,  they  undergo  an  expansion  process. 

The  lumped  parameter  model  developed  by  Coffee  does  permit  a  two-phase 
mixture  and  droplet  combustion  in  the  barrel  region,  i.e.,  assumption  2  is 
relaxed,  but  this  option  has  not  been  exercised  extensively.  The  fully 
one -dimensional  IB  model  recently  developed  by  Gough5  treats  the  combustion 
chamber  as  a  non-homogeneous  region  and  also  permits  a  two-phase  mixture  and 
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droplet  combustion  in  the  barrel  region.  Investigations  using  this  model  of  the 
effects  of  relaxing  assumptions  1  and  2  are  underway,  but  results  are  not  yet 
available . 

Following  the  completion  of  piston  motion,  the  remaining  propellant  in  the 
combustion  chamber  (and  barrel)  is  quickly  consumed,  leading  to  the  "all-burnt" 
condition.  Prior  to  all-burnt,  the  gases  required  to  maintain  the  operating 
pressure,  under  the  conditions  of  rapid  expansion  resulting  from  piston  motion 
in  the  chamber  and  projectile  motion  in  the  barrel,  are  supplied  by  propellant 
combustion,  primarily  in  the  chamber  region.  Upon  burnout,  the  combustion  chamber 
pressure  rapidly  decreases  and  a  rarefaction  wave  would  be  expected  to  move  along 
the  barrel  toward  the  projectile  base. 

This  phenomena  was  originally  suggested  by  Morrison  et  al,6  based  on 
simulations  using  the  one -dimensional  model  developed  by  Gough;2  see  Figure  3. 
(It  appears  similar  to  the  phenomena  reported  in  solid  propellant  systems  utilizing 
stick  propellant  charges  by  Robbins  and  Horst.7)  Recent  analyses  of  data  from 
25 -mm  regenerative  test  firings  conducted  in  the  mid- 1970 's  provide  experimental 
verification  of  the  existence  of  such  a  rarefaction  wave;8  see  Figure  4.  In 
these  tests,  pressure  gages  were  located  at  several  points  along  the  barrel. 
The  resulting  pressure  profiles  are  non-monotonic  along  the  barrel,  with  the 
point  of  maximum  pressure  moving  toward  the  base  of  the  projectile,  indicating 
the  presence  of  a  rarefaction  wave. 

The  ballistic  process  is  completed  with  muzzle  exit  and  "blow-down”  of  the 
gases  remaining  in  the  barrel. 
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IV.  MODIFIED  LAGRANGE  PRESSURE  GRADIENT 


From  the  preceding  discussion,  it  is  obvious  that  the  standard  Lagrange 
pressure  gradient  does  not  accurately  reflect  the  physical  process  in  the  barrel 
region  of  a  regenerative  gun.  The  shortcomings  of  the  standard  Lagrange  model 
are  as  follows: 

1.  The  expansion  of  the  combustion  gases  from  the  chamber  into  the  barrel 
with  a  large  area  reduction  is  not  treated. 

2.  The  non-zero  gas  velocity  at  the  entrance  to  the  tube  and  the  time 
dependence  of  the  mass  of  combustion  gases  in  the  barrel  are  not  accounted  for. 

3.  The  rarefaction  wave  in  the  barrel  after  all-burnt  is  not  simulated. 

In  order  to  simulate  these  processes  in  the  Interior  ballistics  model, 
modifications  to  the  standard  Lagrange  pressure  gradient  model  are  required. 
In  the  following  development,  it  is  assumed  for  simplicity  that  only  combustion 
gases  enter  the  barrel.  However,  the  development  applies  equally  to  a  two-phase 
mixture  of  combustion  gases  and  propellant  droplets  if  the  assumption  of  uniform 
mixture  density  is  made. 

1.  EXPANSION  OF  COMBUSTION  GASES  FROM  THE  COMBUSTION  CHAMBER  INTO  THE  BARREL 


The  barrel  region  is  treated  as  a  single  control  volume  into  which  combustion 
gases  flow  from  the  chamber.  The  left-hand  boundary  of  the  control  volume  is 
defined  by  the  barrel  entrance,  while  the  right  boundary  is  defined  by  the  base 
of  the  projectile.  The  subscript  c  denotes  combustion  chamber  conditions,  L 
denotes  the  left  boundary  conditions  and  R  denotes  the  right  boundary  conditions. 
Therefore , 
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■^ENTRANCE  *  X l  ”  ®  • 


(18) 


*basc  *  x  t "  y  • 


(19) 


and  so  on. 

In  a  steady  state,  isentropic  flow  system,  h  +  v2/ 2  is  a  conserved 
quantity,  where  h  is  the  enthalpy.  Therefore, 


,  1  2  .  1  2  1|7M  2 

ft«*5«.-h|*2Wi  +  2LUr1J  Ui‘ 


(20) 


where  the  last  term  on  the  right -hand- side  has  been  introduced  to  account  for 
the  head  loss  in  the  flow  as  the  gases  enter  the  barrel.  We  note  that 


h  “  cpT  *■  i\P 


(21) 


for  a  gas  obeying  the  Nobel-Abel  equation-of -state  and  assume  that  the  expansion 
of  the  combustion  gases  into  the  tube  is  isentropic,  i.e.. 


'■"’(H-'KiH- 


(22) 


Using  Equations  (21)  and  (22)  in  Equation  (20),  we  obtain,  after  some 
algebraic  manipulation, 
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(23) 


It  is  assumed  ChaC  vc-  0  in  Equation  (23),  which  then  defines  the  pressure  drop 

from  the  chamber  to  the  barrel.  This  equation  is  used  in  both  the  Coffee1-9  and 
Gough2  models  to  connect  the  chamber  and  barrel  regions. 

2.  MODIFIED  LAGRANGE  PRESSURE  GRADIENT;  BEFORE  ALL-BURNT 

We  must  now  solve  the  equations  of  motion  for  the  gases  in  the  barrel, 
Equations  (1)  and  (2),  with  the  boundary  conditions 


x 


t  “  0 . 


(24) 


*„-y. 


(25) 


vL  *  0 . 


(26) 


(27) 


and  with  the  Lagrange  assumption,  Equation  (7).  From  the  continuity  equation, 
we  obtain 
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(8) 


dv  1  dp 
dx  p  dt  ' 


Noting  that  the  mass  of  combustion  gases  in  the  barrel,  m,  is  a  function  of  time 
and  that  p  •m/Alx,  ,  we  obtain 


dp  m  mu  p 
dt  A,x„  A,x2,' 


such  that 


dv  1  <Jp  m  +  up 
dx  p  dt  m  x,  ’ 


(28) 


where  the  time  rate  of  change  of  the  mass  of  the  gases  in  the  barrel,  m,  is  pA,vt. 
Therefore, 


dv  Up  VL 
dx  x,  x,' 


(29) 


Integrating  on  [0,x],  we  have 


and 


(30) 
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(31) 


Substituting  Equation  (30)  and  (31)  in  Equation  (2),  and  noting  that 


(.P  BASE  ^  REs)  w  ' 


Aj 

'  M , 


we  obtain 


+ 


“p-»i 


** 


> 


Integrating  on  [0,x],  we  have 


P(.x)-Pt- p{ 


A, 

2MP 


(  P  BASE 


r  1  /  ^ir 

**  JL 


Note  that  Equations  (30)  and  (33)  reduce  to  (10)  and  (12)  if  vL  -  vL- 

a  -  C. 


Using  the  definition  of  the  space  mean  pressure,  Equation  (14), 
integrate  Equation  (33)  on  [0,xR]  to  obtain 


(32) 


(33) 


0  and 


we  now 
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P’PL- 


m 

nr 


-(p  BAS t  ~  ^Rts) 


m 

3 T, 


(34) 


Similarly,  using  the  definition  for  the  kinetic  energy  of  the  gas,  Equation  (16), 
and  substituting  from  Equation  (30) ,  we  have 


KECas  “-g-Up-Up^-u?)-  (35) 

3.  MODIFIED  LAGRANGE  PRESSURE  GRADIENT;  AFTER  ALL-BURNT 

In  the  description  of  the  RLPG  interior  ballistic  process,  we  discussed  the 
rarefaction  wave  which  travels  along  the  barrel  from  the  chamber  toward  the 
projectile  base  after  propellant  burnout.  We  model  the  rarefaction  wave  as  a 
discontinuity  in  the  spatial  derivatives  of  velocity  and  pressure,  which  travels 
along  the  barrel  at  a  velocity  equal  to  the  average  speed  of  sound  in  the  gas 
plus  the  local  gas  velocity;  see  Figure  5.  Since  this  discontinuity  resembles 
a  "hinge"  in  the  velocity  and  pressure  distributions,  we  have  labeled  it  the 
"hinge  point"  and  quantities  at  the  hinge  point  are  denoted  by  the  subscript  H. 
This  representation  of  a  contact  discontinuity  in  a  gun  tube  was  recently  suggested 
by  Gough.10  The  velocity  of  the  hinge  point,  xB,  is  then 


dt 


vH  +  a. 


(36) 


The  hinge  point  divides  the  barrel  into  two  regions: 
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xR=y 


VR  =  uPROJ 


Figure  5.  Schematic  Showing  the  Hinge-Points  In  the  Bilinear 

Velocity  Distribution  and  Pressure 
Distribution.  Which.,  is  Quadratic  on 
Either  Side .of  the  Hinge-Point. 


Region  I:  0  <  x  5  xH 


and 


Region  II:  xw<x<x,. 


a.  Solution  in  Region  I ;  We  assume  that  the  velocity  distribution  is 
linear  in  both  regions  I  and  II,  i.e.  we  use  the  second  form  of  the  "Lagrange 
Approximation"  noted  by  Comer3.  From  Figure  5,  we  see  that 
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AVl  (x)  _  Vh-Vi 

dx  XH 


0<x<xH, 


(37) 


and  integrating  on  [0,  xH]  we  obtain 


u,(x)-vH 


(38) 


Using  (38)  in  Equation  (2),  we  have 


UP, 
p  dx 


Vl(vh-Vl) 

Xh 


(39) 


We  have  assumed  that  mass  entering  the  barrel  after  all -burnt  is  distributed 
uniformly  over  the  barrel  region,  such  that  there  is  a  uniform  gas  density  in 
the  barrel,  consistent  with  the  basic  Lagrange  assumption.  Integrating  (39)  on 
[0,x]  gives  us  the  pressure  distribution  in  Region  I, 


P,(x)-PL 


where 


(40) 


(41) 
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b.  Solution  in  Region  II :  It  is  assumed  that  at  burnout  xH-xL~0  and 

that  the  hinge  point,  xH,  then  propagates  down  the  barrel  toward  the 
projectile  base  at  a  velocity  given  by  Equation  (36).  Since  xH  propagates  at 
the  speed  of  sound  plus  the  local  gas  velocity,  Region  II  is  acoustically 
isolated  from  events  in  the  combustion  chamber.  Therefore,  we  assume  that  the 
linear  velocity  gradient  in  Region  II  is  determined  by  the  value  of  vL  at  the 
time  of  burnout,  v^,  ,  such  that 


*Vll(.X)  _  Vf-Vu, 

dx  x. 


x,v  £  x  £  x  t . 


(42) 


Then 


and 


(43) 


(44) 


a) 


XH 

x  * 


(45) 


We  now  use  (43)  in  Equation  (2)  to  obtain 


'-‘-r1  ■*.(-) 

P  *x  \x,J 


(46) 
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and  integrating  on  [xa,x],  we  have 


+ 


+ 


( *-*#)“ 


(47) 


where , 


(48) 


(  P  BASE  ~  P  RCS  )  -4  B 

Mo 


(49) 


and  from  (40) ,  we  have 


PH~  PL-pl[("l+VH)Y*hVH~VL)(-VL~a‘>}' 


(50) 


c.  Spare  Mean  Pressure:  After  All -Burnt:  We  again  use  the  definition  of 
Space  Mean  Pressure,  Equation  (14), 


P-  —  f  ‘  P(x)dx -■—([  P,(x)dx*f  P„(x)dx\  . 
XgJo  x i  \yo  J xH  ) 


(51) 


Defining  P;andP,, , 
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-  i  r  * 

P ,  +  —  P /(x)dx . 

X I  J  0 


(52) 


1 

P/i*—  P„(x)dx. 

XtJxH 


we  have 


P, 


g(  vh-vl) 
xh 


uAvh-Vl) 

xh 


\ 

r 


and 


P a  m  P “  "j  '(X X (x *  +  2-x h) 


? 


*  2  vlb  * 


Vg.t.Vt-Vu') 
X < 


(*«-*«)  >• 


Combining  (54)  and  (55)  and  using  Equation  (50),  we  have  for  the  space 
pressure 


P** 


"((-HK-S)] 


V/t 


4- 


X, 


>• 


(53) 


(54) 


(55) 


mean 


(56) 
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Kinetic  Energy  of  the  gas  is  given  by  Equation  (16), 


KE 


CAS 


Substituting  from  (38)  and  (43),  we  have 


KEC 


'mkA,P{fo  v'Cx)dx*fx  v^Cx)dxj- 


(57) 


Defining  KEj  and  KEn ,  we  have 


KE, 


{u2h*vhvl*v2l), 


(58) 


KE„ 


pAlXg 

6 


(59) 


and  finally, 


(60) 
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V.  INTERIOR  BALLISTICS  MODEL 


The  modified  Lagrange  gradient  has  been  implemented  in  the  lumped  parameter 
model  developed  by  Coffee.9  This  model  includes  descriptions  of  the  propellant 
reservoir,  piston  motion,  propellant  injection,  the  combustion  chamber  including 
propellant  combustion,  mass  transfer  to  the  barrel  region,  gas  flow  in  the  barrel, 
and  projectile  motion.  The  changes  required  to  implement  the  modified  Lagrange 
gradient  involved  only  the  portions  of  the  model  dealing  with  mass  transfer  from 
the  chamber  to  the  barrel  and  gas  flow  in  the  barrel. 

1.  STANDARD  LAGRANGE  GRADIENT  MODEL 


In  the  original  version  of  the  Coffee  model,9  a  standard  Lagrange  gradient 
is  used,  and  the  gas  velocity  at  the  entrance  to  the  barrel  is  ignored.  The 
space  mean  pressure  is  calculated  from  the  lumped  parameter  energy  equation  for 
the  barrel.  Phase  is  determined  from  equation  (15) ,  substituting  the  instantaneous 
mass  of  the  gas  in  the  barrel,  m,  for  the  charge,  C, 


^*8ASE 


(61) 


Similarly,  PL  (which  is  equivalent  to  />B«EECHin  this  case)  is  obtained  from  Equation 
(13), 


Pl  ~  P  E 


m 

2AT 


(P  BASE  -  ^Ktj)  ' 


The  mass  flux  into  the  barrel  is 


m  -  pLA,vL, 


(62) 


(63) 
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where  vL  is  obtained  from  Equation  (23)  with  vcm0. 


2  n(Pc-PL) 


(64) 


(Since  we  have  neglected  head  loss  in  performing  the  calculations  discussed  in 
subsequent  sections,  we  have  set  v  “  1  in  writing  Equation  (64).)  The  gas 
density  at  the  left  boundary,  p£,  is  not  equal  to  the  mean  gas  density  in  the 
barrel,  but  is  rather  obtained  from  Equation  (22).  The  kinetic  energy  of  the 
gas  in  the  barrel  is  given  by  Equation  (17). 

2.  MODIFIED  LAGRANGE  GRADIENT  MODEL 


In  the  implementation  of  the  modified  Lagrange  gradient,  ordinary  differential 
equations  for  vt  (i.e.  ,  the  ODE's  for  iiL  are  different  before  and  after  all-burnt) 
must  be  introduced  such  that  vL  is  calculated  directly  in  the  solution  of  the 
coupled  ODE  model  equations.  The  pressure  at  the  left  boundary,  P L,  is  obtained 
by  iteratively  solving  the  non-linear  Equation  (23)  using  the  Newton-Raphson 
technique.  As  in  the  case  of  the  standard  Lagrange  gradient  model,  the  space 
mean  pressure  is  calculated  from  the  lumped  parameter  energy  equation  for  the 
barrel  and  the  mass  flux  into  the  barrel  is  defined  by  Equation  (63). 


a.  Model  Equations  Prior  to  All -Burnt:  Prior  to  all -burnt,  the  pressure 
at  the  projectile  base,  /*BA st,  is  determined  using  Equations  (33)  and  (34).  We 
first  evaluate  Equation  (33)  at  x„  to  obtain 


BASE 


_  m  _  .  m  f.  nl(Up- y£) 

“  Pl  2M~p(P*ASl~PKts)~2T,lVL*  x, 


(6S) 


Using  Equation  (65)  in  Equation  (34)  to  eliminate  the  terras  involving 
up,  vL  and  vL.  we  have 
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(66) 


Alternatively,  we  could  have  used  Equation  (65)  in  Equation  (34)  to  eliminate 
P t,  obtaining 


P  BASE 


m 
3 AT 


-  V  /  ♦ - 

6df|_  x. 


/ 


1  + 


(67) 


(which  reduces  to  the  standard  Lagrange  gradient,  Equation  (61),  for  vl-vl-0.) 

The  required  ordinary  differential  equation  for  i>L  is  then  obtained  by  rewriting 
Equation  (65),  i.e., 


The  kinetic  energy  of  the  gas  in  the  barrel  is  defined  by  Equation  (35). 

b.  Model  Equations  After  All-^urat:  Following  the  completion  of  propellant 
combustion,  the  pressure  at  the  projectile  base,  P 8ASt,  is  given  by  Equation  (47) 
evaluated  at  xK,  i.e., 


P 


BASE  “ 


-r, . 


where  we  have  made  use  of  Equations  (48)  and  (50)  and  where 


(69) 
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(70) 


Tx 


* 


PXH 

2 


xH 


px , 

2 


Xt  Jl  x,J 


Using  Equation  (56)  to  eliminate  the  term  involving  i>,,  we  obtain,  after  some 
algebraic  manipulation, 


3M  p\  2x ,  2x\y 


where , 


(71) 


(72) 


(We  note  that  for  t/16-i)l6-0andxw-0  Equation  (71)  reduces  to  Equation  (61),  the 

standard  Lagrange  gradient.)  The  ordinary  differential  equation  for  vL  is  then 
obtained  by  rewriting  Equation  (69),  i.e., 


(73) 


The  kinetic  energy  of  the  gas  in  the  barrel  is  defined  by  Equation  (60) . 
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The  resulting  sets  of  coupled  ODE  and  algebraic  equations  are  then  solved 
using  the  implicit  ODE  solver,  EPISODE. 

VI.  DESCRIPTION  OF  TEST  CASE 


In  order  to  assess  the  modified  Lagrange  gradient  model,  a  series  of  test 
cases  was  developed.  Since  the  largest  effect  would  be  expected  in  high  velocity, 
high  charge -to -mass  ratio  situations,  a  hypothetical  120-mm  regenerative  LP  tank 
cannon  used  in  an  earlier  study11  was  selected.  The  charge  mass  and  other  gun 
parameters,  with  the  exception  of  liquid  injection  area,  were  held  constant. 
Three  projectile  masses  were  used  to  provide  a  range  of  velocities,  and  the 
liquid  injection  area  was  adjusted  for  each  projectile  mass  to  produce  a  maximum 
chamber  pressure  of  500  MPa.  A  hypothetical  liquid  gun  propellant  was  also  used 
in  the  simulation.  This  propellant  has  the  thermochemistry  of  JA-2  and  the 
physical  properties  of  a  HAN-based  LP.  This  choice  was  made  to  provide  somewhat 
higher  projectile  velocities  than  would  be  achievable  with  a  HAN-based  LP  in  the 
specific  cannon  configuration  used  in  the  study,  in  order  to  permit  evaluation 
of  the  modified  Lagrange  gradient  model  under  the  most  adverse  conditions 
practical.  The  gun  parameters  and  thermochemical  data  used  in  the  simulations 
are  provided  in  Tables  1  and  2 . 


TABLE  1.  120-mm  Cannon  Parameters 


Caliber 

Propellant  Mass 
Projectile  Mass* 
Projectile  Travel 
Initial  Chamber  Volume 
Reservoir  Volume 
Chamber  Diameter 
Liquid  Injection  Area* 
Max  Chamber  Pressure 
Max  Reservoir  Pressure 


120  mm 
16.8  kg 

5  kg,  7  kg,  13  kg 
6.3  m 
6000  cm3 
11700  cm3 
34.2  cm 

94  cm2,  76  cm2,  60  cm2 
500  MPa 
700  MPa 


♦Injection  Area  for  5  kg  projectile  is  94  cm2,  etc. 
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TABLE  2.  Propellant  Data 


Liquid  JA-2  LGP  1846 


Impetus  1140  J/g  898  J/g 

Flame  Temperature  3409  K  2468  K 

Ratio  of  Specific  Heats  1.225  1.223 

Co-Volume  .996  cm3/g  .677  cm3/g 

cp  1.821  J/g-K  1.999  J/g-K 

Density  1.43  g/cm3  1.43  g/cm3 


VII .  RESULTS 

Computer  simulations  for  the  three  test  cases  were  made  using  the  lumped 
parameter  IB  model  developed  by  Coffee9  with  both  the  standard  and  modified 
Lagrange  gradient  models.  Simulations  were  also  run  using  an  IB  model  developed 
by  Gough.2  The  Gough  model  treats  the  reservoir  and  chamber  as  lumped  parameter 
regions  and  utilizes  Equation  (23)  to  connect  the  chamber  with  the  barrel,  which 
is  modeled  as  a  one -dimensional  flow  region. 

Since  a  one -dimensional  model  should  provide  a  more  accurate  simulation  of 
the  barrel  region,  the  Gough  model  was  selected  as  a  baseline  for  evaluation  of 
the  Lagrange  gradient  models . 

Both  models  were  run  assuming  instantaneous  burning  of  propellant  injected 
into  the  combustion  chamber.  As  a  result,  the  simulations  of  the  reservoir  and 
chamber  regions  are  nearly  identical  in  the  Gough  and  Coffee  models ,  and  any 
differences  in  the  overall  model  calculations  are  primarily  due  to  differences 
in  the  simulation  of  the  barrel  region. 

The  calculated  muzzle  velocities  for  the  nine  cases  simulated  (3  projectile 
masses  and  3  gradient  models)  are  presented  in  Table  3. 
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TABLE  3. 


Comparison  of  Calculated  Muzzle  Velocities 


Gough  Model 
One - Dimens ional 
Flow 


Coffee  Model 
Standard  Lagrange 
Gradient 


Coffee  Model 
Modified  Lagrange 
Gradient 


jj  ectile 

Muzzle 

Muzzle 

Muzzle 

Mass 

Velocity 

Velocity 

Difference 

Velocity 

Difference 

(kg) 

(m/s) 

(m/s) 

(») 

(m/s) 

(%) 

5.0 

2020 

2121 

5.0 

1952 

-3.4 

7.0 

1894 

1983 

4.7 

1849 

-2.3 

13.0 

1563 

1622 

3.8 

1546 

-1.1 

The  differences  between  the  muzzle  velocities  calculated  using  the  Lagrange 
gradients  and  those  obtained  from  the  Gough  models  are  presented  as  a  percent 
of  the  baseline  velocity.  As  can  be  seen,  use  of  the  standard  Lagrange  gradient 
consistently  results  in  muzzle  velocities  higher  than  the  baseline  case,  while 
use  of  the  modified  Lagrange  gradient  results  in  velocities  which  are  lower  than 
the  baseline.  Overall,  the  velocities  obtained  using  the  modified  Lagrange 
gradient  are  somewhat  closer  to  the  baseline  than  those  obtained  with  the  standard 
gradient,  and  the  agreement  degrades  with  increasing  velocity. 


The  pressure  versus  time  curves  from  the  simulations  with  the  7  kg  projectile 
mass  are  presented  in  Figures  6a  and  6b.  In  Figure  6a,  the  pressure  versus  time 
curves  at  three  locations  (chamber,  barrel  entrance  and  projectile  base)  obtained 
using  the  Coffee  model  with  the  standard  Lagrange  gradient  (dotted  line)  are 
presented  along  with  the  corresponding  pressure  curves  obtained  using  the  Gough 
model  for  comparison.  A  similar  comparison  is  presented  in  Figure  6b  for  the 
modified  Lagrange  gradient. 

As  can  be  seen  in  Figure  6a,  the  chamber  pressures  are  quite  similar,  but 
the  pressure  versus  time  curves  at  the  barrel  entrance  and  projectile  base  show 
significant  differences.  The  standard  Lagrange  gradient  results  in  a  pressure 
at  the  projectile  base  which  is  consistently  higher  than  that  obtained  with  the 
one -dimensional  model,  resulting  in  a  higher  muzzle  velocity;  see  Table  3.  The 
pressure  at  the  barrel  entrance  for  the  standard  Lagrange  gradient  case  is 
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In  comparison,  the  modified  Lagrange  gradient  model  results  show  very  good 
agreement  with  the  one -dimensional  model;  see  Figure  6b.  The  chamber  and  barrel 
entrance  pressures  are  almost  indistinguishable  over  the  entire  ballistic  cycle 
for  the  two  simulations.  The  projectile  base  pressures  are  also  nearly 
indistinguishable  over  the  first  half  of  the  ballistic  cycle,  but  the  modified 
Lagrange  curve  Increasingly  departs  from  the  baseline  curve  after  about  4.0  ms. 
This  departure  is  below  the  baseline,  consistent  with  the  velocity  difference 
noted  in  Table  3.  This  discrepancy  between  the  modified  Lagrange  and 
one -dimensional  results  is  attributed  to  the  decreasing  validity  of  both  the 
basic  Lagrange  assumption,  uniform  gas  density,  as  the  projectile  becomes  farther 
removed  from  the  barrel  entrance,  and  the  assumption  of  fixed  vLb  and  vLb  after 
all -burnt. 

Pressure  and  velocity  profiles  in  the  barrel  region  from  the  simulations 
for  the  7  kg  projectile  using  the  modified  Lagrange  and  one -dimensional  models 
are  presented  in  Figures  7-11.  The  times  for  which  these  profiles  were  plotted 
were  chosen  to  coincide  with  maximum  pressure  at  the  projectile  base  (3.3  ms), 
maximum  chamber  pressure  (4.7  m/s),  a  time  just  after  all-burnt  (6.1  ms),  a  time 
when  the  rarefaction  wave  has  propagated  some  distance  down  the  barrel  (6.5  ms), 
and  muzzle  exit  (7.4  m/s). 

Prior  to  all-burnt  (3.3  ms  and  4.7  ms),  the  modified  Lagrange  model  produces 
velocity  and  pressure  profiles  which  have  the  correct  shape  and  are  quite  close 
to  the  one -dimensional  profiles  over  the  entire  barrel  region;  see  Figures  7  and 
8.  The  departure  in  pressure  at  the  projectile  base,  noted  earlier,  as  well  as 
a  corresponding  difference  in  velocity  near  the  projectile  base,  has  begun  to 
appear  in  Figure  8 . 
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In  Figure  9,  Che  rarefaction  wave  can  be  seen  just  beginning  to  propagate 
along  the  barrel  toward  the  projectile  base.  Again,  the  overall  shape  and 
magnitude  of  the  velocity  and  pressure  profiles  are  in  good  agreement.  The 
differences  in  velocity  and  pressure  near  the  projectile  base  have  become  more 
pronounced,  about  6-7%  in  pressure  at  6.1  ms  versus  about  4%  at  4.7  ms. 


.1  0.5 


Figure  8 . 


Slightly  later  in  time,  at  6.5  ms,  the  rarefaction  wave  is  well  developed 
and  has  propagated  about  one  quarter  of  the  distance  to  the  projectile  base;  see 
Figure  10.  The  velocity  and  pressure  profiles  are  still  in  good  agreement  with 
the  baseline  case,  and  the  hinge  point  is  in  approximately  the  correct  position. 


33 


Figure  9.  Comparison  of  Velocity  and  Pressure  Distribution  in  the 

Barrel  at  ,t  -  6,1  as.  from  Simulations 
Using  Gough  Model2  (Solid  Line>  and 
the  Modified  Lagrange  Gradient  (Dashed  Line) 
for  the  7  kg  Projectile. 


The  pressure  and  velocity  differences  near  the  projectile  base  have  continued 
to  increase;  the  pressure  difference  at  the  projectile  base  has  increased  to 
about  8% . 


At  muzzle  exit,  Figure  11,  the  velocity  profiles  are  still  in  very  good 
agreement,  but  the  modified  Lagrange  pressure  profile  has  deteriorated  somewhat. 


3' 


0.0  0.5  1.0  1.5  Z.o  2.5  3.0  3.5  4.0  4.5  5.0 


x  (m) 


0.0  0.5  1.0  1.5  2.0  2.5  J.O  J.5  4.0  4.5  5.0 


x  (m) 


Figure  10. 


The  shape  of  the  pressure  profile  is  still  approximately  correct;  however,  the 
departure  from  the  baseline  profile  is  apparent,  the  pressure  difference  at  the 
projectile  base  having  increased  to  approximately  10%.  This  discrepancy  does 
not  reflect  the  overall  accuracy  of  the  approximation,  however,  since  the 
cumulative  difference  in  projectile  velocity  is  only  about  2.5%  in  this  case. 


Figure  11 .  Comparison  of  Velocity  and  Pressure  Distribution  in  the 

Barrel  at  t  -  7.4  ms.  from  Simulations 
Using  Gough  Model?  (Solid  Llnel  and 
the  Modified  Lagrange  Gradient  (Dashed  Line1) 
for  the  7  kg  Proj-ectlla. 


VIII.  CONCLUSIONS 


A  modification  of  the  standard  Lagrange  gradient  model  has  been  developed 
to  simulate  the  gas  flow  in  the  barrel  of  a  regenerative  liquid  propellant  gun. 
This  model  accounts  for  the  non-zero  gas  velocity  at  the  entrance  to  the  barrel, 
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and  the  existence  of  a  rarefaction  wave  which  travels  along  the  barrel  from  the 
chamber  toward  the  projectile  base  after  all -burnt.  The  modified  Lagrange 
gradient  has  been  implemented  in  the  lumped  parameter  interior  ballistic  model 
developed  by  Coffee,  and  comparisons  have  been  made  among  the  resulting  IB  model, 
the  Coffee  model  with  a  standard  Lagrange  gradient,  and  the  Gough  model  with  a 
one -dimensional  barrel  flow  model  Computer  simulations  of  a  high  velocity 
120-mm  tank  cannon  have  been  run,  using  these  three  models,  for  3  projectile 
masses . 

The  modified  Lagrange  gradient  model  resulted  in  muzzle  velocities  1.1-3. 4% 
lower  than  the  one -dimensional  model  over  the  velocity  range  1563-2020  m/s,  while 
the  standard  Lagrange  model  produced  muzzle  velocities  3. 5-4. 7%  higher  than  the 
one -dimensional  model.  The  pressure  versus  time  curves  obtained  using  the 
modified  Lagrange  gradient  are  in  very  good  agreement  with  those  from  the 
one -dimensional  model,  although  some  discrepancy  is  observed  at  the  projectile 
base  over  the  latter  half  of  the  ballistic  cycle.  A  similar  comparison  for  the 
standard  Lagrange  gradient  case  shows  very  poor  agreement  with  the  one-dimensional 
model  at  the  barrel  entrance  and  at  the  projectile  base  over  most  of  the  ballistic 
cycle . 

Detailed  comparisons  of  the  velocity  and  pressure  profiles  in  the  barrel 
from  the  modified  Lagrange  and  one -dimensional  models  at  five  different  times 
during  the  ballistic  cycle  have  been  made.  These  comparisons  show  that  the 
profiles  obtained  using  the  modified  Lagrange  gradient  have  the  correct  shape 
and,  with  the  exception  of  the  pressure  profile  near  muzzle  exit,  are  in  excellent 
agreement  with  the  one -dimensional  results.  Overall,  the  modified  Lagrange 
gradient  model  appears  to  be  capable  of  accurately  representing  the  physical 
processes  in  the  barrel  region  of  a  regenerative  liquid  propellant  gun  and  of 
providing  an  excellent  overall  simulation  of  the  IB  process,  even  in  the  high 
velocity  regime. 
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Subscripts : 


BASE 

BREECH 

c 

H 

L 

Lb 

R 


Projectile  Base 

Gun  Breech 

Combustion  Chamber 

Hinge-Point 

Barrel  Entrance 

Barrel  Entrance  at  All-Burnt 

Projectile  Base 
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